Statistical Reconstruction of Qutrits 
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We discuss a procedure of measurement followed by the reproduction of the quantum state of a 
three-level optical system - a frequency- and spatially degenerate two-photon field. The method of 
statistical estimation of the quantum state based on solving the likelihood equation and analyzing the 
statistical properties of the obtained estimates is developed. Using the root approach of estimating 
quantum states, the initial two-photon state vector is reproduced from the measured fourth moments 
in the field . The developed approach applied to quantum states reconstruction is based on the 
amplitudes of mutually complementary processes. Classical algorithm of statistical estimation based 
on the Fisher information matrix is generalized to the case of quantum systems obeying Bohr's 
complementarity principle. It has been experimentally proved that biphoton-qutrit states can be 
reconstructed with the fidelity of 0.995-0.999 and higher. 
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I. INTRODUCTION 

The ability of measuring quantum states is of funda- 
mental interest because it provides a powerful tool for the 
analysis of basic concepts of quantum theory, such as the 
fundamentally statistical nature of its predictions, the su- 
perposition principle, Bohr's complementarity principle, 
etc. To measure quantum state one needs to perform 
some projective measurements on the state and then to 
apply some computation procedure to the data. The first 
step is a genuine measurement consisting of a set of op- 
erations on the representatives of a quantum statistical 
(pure or mixed) ensemble. As a result of such operation 
an experimentalist acquires a set of frequencies at which 
particular events occur. In the second step a mathemati- 
cal procedure is applied to the statistical data obtained in 
the previous step to reconstruct the quantum state. The 
present paper is devoted to the state reconstruction for 
the optical three-level systems. The object under study 
is the polarization state of a frequency- and spatially de- 
generate biphoton field Q . 

The necessity of the adequate measurement of the 
states of such systems is caused not only by fundamental 
interest but also by some applications. For example, it 
has been shown that the security of the key distribution 
in quantum cryptography is associated with the dimen- 
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sionality of the Hilbert space for the states in use 0]. 
From this point of view certain hopes are pinned on the 
three- level systems or qutrits |^^|^ rather than qubits. 

We should mention that there are other implementa- 
tions of three-level optical systems. The most familiar 
ones deal with three-arm interferometers and lower- 
order transverse spatial modes of optical field, realized 
with holograms 0, 0, 0|- Polarization-entangled four- 
photon fields, which are equivalent to two entangled spin- 
1 particles were studied in [lo| . 

The paper is organized as follows. In Sec. II we discuss 
the main properties of qutrits based on the polarization 
state of biphoton field. We focus on their preparation, 
visual representation on Poincare sphere, unitary trans- 
formation by phase plates. Then we consider the coher- 
ence matrix, which characterizes completely the prop- 
erties of biphotons-qutrits in the fourth field moments. 
Sec. Ill is devoted to the methods of biphotons-qutrits 
measurement, in particular, we introduce two quantum 
tomography protocols and discuss in detail their exper- 
imental implementation. We conclude this part with 
the analysis of statistical reconstruction for qutrits from 
the outcomes of mutually complementary measurements. 
Sec. IV deals with the methods of quantum state re- 
construction. Namely we consider the least-squares and 
maximum-likelihood methods and apply these tools to 
analysis the data obtained in quantum tomography. In 
Appendix we explore the problem of statistical fluctua- 
tions of the state vector which is important for the esti- 
mation and control of precision and stability of quantum 
information. 



2 



II. QUTRITS BASED ON BIPHOTONS 

A. Preparation 

Biphoton field is a coherent mixture of two-photon 
Fock states and the vacuum state flU : 



ICjl 



^ = \vac) 
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kski 



(1) 



where , l^r ) denotes the state with one (signal) pho- 
ton in the mode ks and one (idler) photon in the mode 
ki . The coefficient F^- ^r. is called the biphoton amplitude 
[T^ . because its squared modulus gives a probability to 
register two photons in modes ks and ki. 

Let us consider the collinear and frequency degenerate 
regime, for which k^ ^ ki, uig ~ uii and ujs + uji = ujp, 
where tOp is the laser pump frequency. We further restrict 
our discussion to biphotons that are indistinguishable in 
terms of spatial, spectral, or temporal parameters. From 
the point of view of polarization there are three natural 
states of biphotons, namely, ^'i = |2,0), ^2 = 
and 1-3 = |0,2). Here the notation |2,0) = |2jj,0y), 
for example, indicates that there are two photons in the 
horizontal (H) polarization mode, while no photons are 
present in the orthogonal vertical (V) mode. These basic 
states can be generated using type-I (for 'I'l and "if 3 ) 
and type-II (for ^'2) phase matching. Since only two- 
photon Fock states arc considered, for the state |m,n) 
the condition m + n — 2 must be satisfied. 

Any arbitrary pure polarization state of biphoton field 
can be expressed in terms of three complex amplitudes 
ci, C2, and C3: 



|c)=ci|2,0)+C2|l,l)+C3|0,2), (2) 



where Cj = |cj| exp{i(/?j}, ^ |cjp = 1. The vector |c) — 

i=i 

(ci, C2, C3) represents a three-state state or qutrit. 

There is an important note concerning the state- vector 
(2). In principle, one can write the complete polarization 
state in the form 



|c) =Ci|2ff,0y) -f C2|lv,lff) +C3 

where the terms and might be distin- 

guishable somehow, for example, if the photon with ver- 
tical polarization comes first with respect to the photon 
with horizontal polarization. However we consider par- 
ticular two-mode polarization state so photons differ in 
polarization only and there are no other parameters re- 
sponsible for their distinguishability. 

In general, to generate an arbitrary qutrit state one 
needs to put three nonlinear crystals separated in space 
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FIG. 1: Preparation of an arbitrary qutrit based on bipho- 
tons, in principle. Three nonlinear crystals placed in the com- 
mon pump generate biphotons with type-I (1, 3) and type-II 
(2) phase-matching. Three attenuators (|ci|^, |cip, |c3p) and 
three phase-shifters {ipi , ip2 , 'fia) allow one to control three 
complex amplitudes Ci, C2 and C3. 



into a common pump and superpose the biphoton fields 
generated by the three crystals coherently or incoherently 
(Fig. I). 



B. Representation of qutrits using the Poincare 
sphere 

Sometimes it is very convenient to use visual represen- 
tation of the state. For example a single-photon pure po- 
larization state (qubit) may be mapped onto the Poincare 
sphere (3-dimensional Euclidian sphere). A (pure) qubit 
state is determined by polar and azimuthal angles (1?, </>) 
in spherical coordinates. Any unitary polarization trans- 
formation of the qubit is represented by the correspond- 
ing rotation of the sphere. Thus, in order to learn the 
final transformed state one just has to apply the rotation 
operation using certain rules. 

It would be helpful to use same visual representation of 
a qutrit using the Poincare sphere. Although the gener- 
alization of the Poincare sphere for qutrits has been dis- 
cussed earlier ^3 suggest an alternative approach, 
which allows us to manipulate with qutrits in natural 3- 
D space rather than in sophisticated 8-D space. Let us 
map the polarization state of a biphoton into a pair of 
points on the sphere (but this is not the two-qubit case 
since the states 1 77, and \V,H) arc indistinguishable). 
In this representation each photon forming the biphoton 
is plotted as a single point on the Poincare sphere, so the 
qutrit state vector is represented by 
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where a^(t?j, 4>i) and at(?9j, 4>s) E^rc the creation operators 
in idler and signal polarization modes and a^(i?m, 4>m) = 
cos(??m/2)a'l' -I- e"^" sin(??m/2)6^, m=i,s. Note that oper- 
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action of such elements on the state (2) is described by 
the matrix flSll: 



G 



L 
R 



L 
R 
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where 



t — cos (5 + i sin 8 cos 2a, 



isin^ sin2Q;, 



(7) 



(8) 



FIG. 2; Representation of a qutrit using the Poincare sphere. 
Figures a, b, c show three basic states forming superposition 
(2). Figure d represents the state of an arbitrary qutrit. 



ators = a|^,6t = ay are creation operators for H— 
and 1/— polarized photons. 

It is well-known that the number of real parameters 
characterizing a quantum state is determined by the di- 
mension of the Hilbert space (s). For a pure state, 



6 = 7r(7io — ne)h/X is the optical thickness of the plate, 
h is its geometrical thickness, a is the orientation angle 
between the optical axis of the plate and one of the basis, 
for example, vertical direction. 

Let us consider the action of the half-lambda plate 
on a particular state — --^(|2,0) — |0,2)), when the 
plate is oriented at 22.5°. For the state there are two 
nonzero amplitudes ci = C3 = 1/V2 and there is only one 
relative phase (^13 = (pi ip^ ~ n. Taking into account 
that for a half-lambda plate S — 7r/2, the corresponding 
transmission and reflection coefficients are 



^pure — 2s 2, 



and for mixed states. 



mixed 



(5a) 



(5b) 



t 

Thus the matrix G has the form 



(9) 



G 




According to (5a, b), four real parameters determine com- 
pletely the pure state of a qutrit, so in the Poincare 

sphere representation these parameters are simply the y 
four spherical angles {-di, (j)i;^s, 'Ps)- The links be- 
tween the angles {'dj, (jn; i^s, (jis) and the amplitudes Cj = Hence, acting by matrix G on the state we get 
\cj I exp itpj are derived in |l4l | . As an example three basic 

states = |2, 0), 'I'2 = |1, 1), and = |0, 2) are shown q 
in Fig. 2. It can be shown that the polarization degree G^ \ = 



(10) 



of a qutrit P = V|cip - |c3p 2|c|c2 -f- c^cgp Q has 
a clear geometrical meaning: it is defined by the angle (3 
between the pair of points on the Poincare sphere as seen 
from its center: 



V2 



1 =*2 
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2cos(/3/2) 
1 + cos2(/3/2)' 



(6) 



For the states '^i and '^^ the polarization degree takes 
values Pi_3 = 1, since two points coincide on the sphere 
and /3 = 0. For the second state, ^"2, two points are 
positioned at the opposite sides of the sphere, that is 
why (3/2 = 7r/2 and P2 = 0. 



C. Transformation 

Experimentally a unitary transformation of the polar- 
ization state (2) can be achieved by placing any retarda- 
tion plates, rotators etc. into the biphoton beam. The 



Note that such kind of transformations cannot change the 
polarization degree of a qutrit. For the state 'I'_l chosen 
above, as well as for the state ^"2, the polarization degree 
P is zero. 

In the experiment described below we used a simpler 
way to generate qutrits. Biphotons were produced via 
collinear frequency-degenerate spontaneous parametric 
down conversion in a nonlinear crystal (BBO, type-I or 
type-II phase matching). For type-I phase matching the 
polarization of both created photons was vertical; i.e., the 
state ^3 was generated. Then, this state was transformed 
using a quartz plate with a fixed optical thickness. By 
changing the angle of the plate, the state ^'3 ~ |0/f, 2y) 
is transformed according to the formula |ci„) = G'^z- 
For the case of type-II phase matching the final state is 
|cm) = G^'2. Of course the state |ci„) does not involve 
all possible qutrit states because the transformation given 
by matrix (7) preserves the polarization degree. Anyway 
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using such a transformation, we select some subset of 
qutrits to work with. 

Such a simple method of the state prepara- 
tion/transformation was chosen in order to be able to 
compare the results of reconstruction with the parame- 
ters of the input states, which should be known with a 
high accuracy. The purpose of this work is the recon- 
struction of the initial state |cm). 

We would like to emphasize that only pure qutrit states 
are accessible by this method. To create a mixed state, 
some more complicated method is to be used. This 
method allows one to create arbitrary qutrit states and 
it implies a possibility to introduce controlled delay be- 
tween three fundamental states forming the qutrit which 
could exceed the coherence length of the laser pump . 

D. Coherence matrix 

We introduced only qualitative description of the 
qutrits based on biphotons so far. The quantitative 
measure characterizing the polarization properties of any 
single-mode state in the forth moment in the field (includ- 
ing biphoton state) was proposed by D.Klyshko in [T7l |. 
It is a matrix consisting of six fourth-order moments of 
the electromagnetic field. An ordered set of such mo- 
ments can be obtained using the direct product of 2 x 2- 
coherence matrixes for both qubits. After normal order- 
ing, averaging, and crossing out the redundant row and 
column the matrix takes the following form: 

( A D E\ 
Ki= I D* C F \ . (11) 
\E* F* B ) 

The diagonal elements are formed by real moments, 
which characterize the intensity correlation in two po- 
larization modes H and V: 

A={a'^'^a^), B={bH^), C = {a^ b'' ab) , (12) 
Nondiagonal moments are complex: 

D={a^^ab), E={aH''), F={a^W). (13) 

Three real moments (12) and three complex ones (13) 
completely determine the state under consideration. The 
elements of the matrix (11) are expressed through the el- 
ements of the polarization density matrix. The normal- 
ization condition, 

A + B + 2C = 2, (14) 

reduces the number of independent real parameters, so 
for a mixed state we get 8 parameters as expected. In the 
special case of a pure biphoton state, taking the average 



in Eqs. (12,13) over the state (2), we obtain the matrix 
components in the following form: 

A = 2\ci\\B = 2\c^\\C = \c2\\ (15) 



D:^V2clc2, £; = 2c*C3, F = V2cJc3 (16) 

So the links between the polarization density matrix and 
the matrix (11) can be found comparing the correspond- 
ing components of {K4)rnk and of p = |c)(c|; pmk = 
CmCk*]'m,k = 1,2,3 for a pure state and pmk = CmCk* 
for a mixed state where the averaging, as usual, is taken 
over the classical probability distribution. The statistics 
of the field is assumed to be stationary and ergodic so 
the time-averaged values of the observed quantities can 
be described in terms of a quantum statistical ensemble. 
In this case (...) = Tr{p...), where p is the polarization 
density operator. 

III. METHODS OF MEASUREMENT 

What does it mean to measure the unknown state (2)? 
From the experimental point of view, it means that the 
experimentalist has to measure a complete set of real pa- 
rameters (moments) determining the state. To do this 
the state must be subject to a set of unitary polariza- 
tion transformations and projective measurements. By 
doing this one picks out the outcomes, which are pro- 
portional to the corresponding moments (12, 13) or their 
linear combination. This procedure is known as quan- 
tum tomography. The quantum state can be represented 
using either the wave function, density matrix, or quasi- 
probability function (Wigner function). Probably the 
correct way to use the term "quantum tomography" is 
only for the reconstruction of the quasi-probability func- 
tion because it gives the graphical representation of the 
state as a 3D plot. Nevertheless the term "quantum to- 
mography" is also used for a general procedure of com- 
plete state reconstruction. For a brief review among the 
papers where this procedure was realized experimentally, 
let us mention the works 0, 0, |23 related to states 
defined by continuous variables. For states character- 
ized by discrete variables, such as two polarization-spatial 
qubits, quantum tomography was realized in [23|. Re- 
cently quantum tomography has been performed for or- 
bital angular momentum entangled qutrits Q etc. 

The physical idea behind the tomography procedure is 
performing measurements of a ppr opriately complete set 
of observables called quorum j22| or just "looking" at 
the state from different positions. The minimal number 
of such positions might be the number of real parameters 
determining the state. 

According to Bohr's complementarity principle, it is 
impossible to measure all moments (12,13) simultane- 
ously, operating with a single qutrit only. So to perform 



a complete set of measurements one needs to generate a 
lot of representatives of a quantum ensemble. 

First of all, let us mention that at present, the only 
realistic way to register single-mode biphoton field is us- 
ing the Brown- Twiss scheme. This scheme consists of a 
beam-splitter followed by a pair of detectors connected 
with the coincidence circuit. It means that registration 
of a single biphoton, which carries the state (2), can give 
only a single event at the output of the experimental set- 
up with some probability. So the statistical treatment 
of the outcomes becomes extremely important. For cor- 
relations between polarization degrees of freedom, which 
is essential in the case under consideration, the Brown- 
Twiss scheme must be accomplished with polarization 
filters introduced into each arm. 



A. Qutrit tomography protocols 

We proposed two methods to perform polarization re- 
construction of a biphoton qutrit state |ci„). 
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FIG. 3: Measurement block for Protocol 1. The Brown- Twiss 
scheme for measuring intensity correlation between two polar- 
ization modes. After spatial separation at the non-polarizing 
beam splitter (BS), signal (s) and idler (i) photons propagate 
through the quarter- and half wave plates, polarizing prisms 
(V), focusing lenses (F) and interference filters (IF) in two 
channels. Finally, photons are registered by detectors (D). 
The coincidence rate from the output of the coincidence cir- 
cuit (CC) is proportional to the fourth moment in the field 

{Rs^)■ 



Protocol 1. 



state vector of the signal (idler) photon; 



The idea of the first method is splitting the state \cin) 
into two spatial modes and performing transformations 
over two photons independently (Fig 3). These transfor- 
mations can be done using polarization filters placed in 
front of detectors. Each filter consists of a sequence of 
quarter- and half-wave plates and a polarization prism, 
which picks out definite linear polarization, for example, 
the vertical one. A narrowband filter centered at the dou- 
bled pump wavelength A = 2Ap serves to make biphotons 
emitted from different sources indistinguishable in fre- 
quency as well as to reduce the background noise. An 
event is considered to be detected, if a pulse appears at 
the output of the coincidence circuit. Approximately in 
half of trials, one of the photons (signal, by convention) 
forming a biphoton is going to one of the detectors, while 
the other one (idler) is going to the other detector. In the 
remaining cases, both photons appear in the same out- 
put beam-splitter arm, and these events are not selected 
because they do not contribute to coincidences. 

In the Hcisenbcrg representation the polarization 
transformation for each beam-splitter output port is 
given by: 









n 




= 


[I 





7r/2,0) 



y.Dx/A{5 = 7r/4,x) 




(17) 



Four 2x2 matrixes in the right-hand side of Eq. (17) 
describe the action of the non-polarizing beam-splitter, 
A/4-, A/2- plates and vertical polarization prism on the 



D 



A/2, A/4 



t r 

-r* t* 



where r and t - are the coefficients introduced in Eq. (8), 
so for a A/4-platc {5 ~ 7i'/4), 



1 i 
= ;^(1 + «cos2x),rA/4 = -^sin2x (18a) 



and for a A/2-plate {6 = 7r/2) 



t\/2 = * cos(26'), r;^/2 = i sin(26'). 



(18b) 



Thus, there are four real parameters (two for each 
channel) that determine polarization transformations. 
Namely, these parameters are orientation angles for two 
pairs of wave plates: 0i,Xii^2,X2- 

As it was mentioned above, the output of the Brown- 
Twiss scheme is the coincidence rate of the pulses coming 
from two detectors Dg and Di. The corresponding mo- 
ment of the fourth order in the field has the following 
structure: 

Rs,, cx {h'J'b'.h'gb',) = i?(0i, XI, ^2, X2) (19) 

In the most general case this moment contains a linear 
combination of six moments (12, 13) forming the ma- 
trix K^. So the main purpose of the quantum tomogra- 
phy procedure is extracting these six moments from the 
set-up outcomes by varying the four parameters of the 
polarization Brown- Twiss scheme. 
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Consider some special examples, which give the corre- 
sponding lines in the complete protocol introduced below 
(Table I). 

First of all. it is obvious that for measuring real mo- 
ments (12) one needs to make polarization filters transmit 
both photon with horizontal polarizations to measure A, 
both photons with vertical polarization to measure B and 
one photon with vertical and another one with horizontal 
polarizations to measure C. To do this all quarter-wave 
plates should be oriented at zero degrees, then to install 
both half- wave plates at zero degrees for measuring B; at 
9s = 45° and 9i = 45° for measuring A; and at 9s = 0°, 
9i = 45° for measuring C. These settings pick out the 
squared modulus of corresponding amplitudes C3 , ci and 

C2- 

The next example shows how to measure one of the 
complex moments (13). To measure the real part of the 
moment D, let us set the wave-plates in the Brown-Twiss 
scheme in the following way. 
The idler channel: 



A/4:x» = 0°,i?A/4 = 



1 

V2 



1 + i 






1-i 



\/2:9,^i5°,D^/2 
The signal channel: 



i 

1 



(20a) 



(20b) 



A/4:x. = 45°,i?A/4 



V2 U 1 



A/2 : 0, =22.5°, Da/2 = 



(21a) 



(21b) 



Substituting these matrices into Eq. (17) and taking into 
account the commutation rules for the creation and an- 
nihilation operators it is easy to get the final moment to 
be measured: 



R = {c\blblbsb,\c) = 1/8{A + C- 2ReD). 

A complete set of the measurements called the tomog- 
raphy protocol is presented in Table I. Each row corre- 
sponds to the setting of the plates to measure the moment 
placed in the sixth column. The last one corresponds to 
the amplitude of the process (see below). 

This protocol was suggested and developed in [23L |23 . 
A similar protocol was considered in details earlier |2l| 
for estimating polarization state of a biphoton field, gen- 
erated in a frequency degenerate non-coUinear mode. In 
this case the biphoton field is represented as a pair of 
polarization qubits. 

Before describing the second method of the state mea- 
surement let us make some remarks. 



*We assume that the source generating qutrirs is sta- 
tionary. Since each measurement eliminates a qutrit one 
has to be sure that there are a lot of copies of the initial 
state; each copy must be prepared in the same quantum 
state. Such ensemble approach guarantees that the ex- 
perimentalist deals with the same quantum state in all 
trials. In other words, the outcomes provide him with 
the information about the same quantum state and elim- 
ination of a particular state does not affect the rest ones. 

**The outcomes of the set-up are numbers related to 
the corresponding moments (19). Usually this number 
is the coincidence counting rate or the number of coinci- 
dences in a fixed time interval. Due to the necessity of a 
proper normalization of the state under investigation, the 
number of independent real parameters grows up. The 
normalization is obtained from the measurement of mo- 
ments A, B and C. Furthermore, only cosine and sine of 
the phases ipi2 and ipis can be measured in experiment 
as there is no way to measure the phases directly. That 
is why the final number of real parameters to be mea- 
sured in experiment is 7 for a pure qutrit state and 9 for 
a mixed state. 

*** To minimize the errors caused by independent sta- 
tistical fluctuations of the outcomes, the number of mo- 
ments (12, 13) entering in Eq. (19) should be minimal. 



2. Protocol 2 

In the second method of quantum tomography, a 
biphoton-qutrit being measured is first subject to a se- 
quence of unitary transformations and, for each of such 
transformation, it is fed to the Brown-Twiss scheme set- 
tled for measuring a fixed moment. Using the wave plate 
with arbitrary optical thickness, one can achieve the quo- 
rum varying the orientation of the plate ^. 

In the most general case the coincidence counting 
rate in this protocol is a periodic function of fj,, more- 
over, its Fourie expansion contains nine harmonics of 
fi : cos(O^), cos(2^), sm{2^), cos(4^), sin(4/.i), cos(6/.i), 
sin(6/i), cos(8/.i), sin(8/^). These harmonics depend lin- 
early on the nine moments A, B, C, ReD, ImD, ReE, 
ImE, KcF, ImF. In other words, there is a 9 x 9 matrix 
T that links these nine harmonics to the nine moments 
as shown below. 



/ COSOyU \ 




( ^ \ 


cos 2^ 




B 


sin 2fj. 




C 


cos 4/i 




YlcF 


sin 4/i 


= T 


ImF 


cos 6fi 




ReZ? 


sin 6/i 




ImD 


cos 8yU 




RcE 


\ sin ByU / 




\ ImD / 



Unfortunately, the inverse matrix does not always exist. 
To simplify the problem we put only a single wave plate 
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TABLE I: Protocol 1. Each line contains orientation of the half (Os^i) and quarter (xs,0 wave plates in the measurement block. 
Last two columns show the corresponding moment 7?„ and the process amplitude Mv{v = 1, ..9). 





Parameters of the experimental set-up 


Moments to be measured 


Amplitude of the process 


V 


Xs 




Xi 








1. 





45° 





-45° 


A/ A 


ci/V2 


2. 





45° 








C/4 


C2/2 


3. 














B/4 


C3/V2 


4. 


45° 











1/8(B + C + 2ImF) 




5. 


45° 


22.5° 








1/8(B + C - 2RcF) 




6. 


45° 


22.5° 





-45° 


1/9.{A + C -2RcD) 




7. 


45° 








-45° 


l/'&iA + C + 2lmD) 




8. 


-45° 


11.25° 


-45° 


11.25° 


l/16(A + B-2Im£) 




9. 


45° 


22.5° 


-45° 


22.5° 


1/I6(y4 + B- 2ReE) 


275=1 ^ 575=^ 



QP2 BS hi V F IF D 




FIG. 4: Measurement block for Protocol 2. Additional control 
quartz plate (QP2) serves as state |ci„) tomography trans- 
former. Only a single wave plate is introduced in each channel. 



with fixed optical thickness {Sg = 7r/4, (5; = tt/2) and 
fixed orientation Xs , Oi in each channel of the Brown- 
Twiss scheme. In order to make sure the inverse ma- 
trix exists one needs to maximize the determinant of the 
matrix T over the orientations of the plates Xs,Oi- Af- 
ter accomplishing this procedure wc obtain Xs ~ 19°, 
0^ « -28.5° (Fig.4). 

Instead of finding the links between the harmonics and 
moments, there is a more elegant method to reconstruct 
the quantum state using the second protocol (see section 
III.D). This method is considered in the present work for 
the first time. 



B. Experimental implementation. Protocol 1 

The experimental set-up for the quantum tomography 
of qutrits using protocol 1 is shown in Fig. 5. The prepara- 
tion block includes a 2-mm BBO crystal with either type- 
I or type-II degenerate and collinear phase-matching, 
which is pumped with cw-argon laser operated at 351nm 
wavelength. In the case of typc-II phase-matching, an ad- 
ditional quartz compensator is introduced right after the 
crystal. The state "i/s = |0,2) (for typc-I) or *2 = |1, 1) 



(for type-II) generated in the crystal is fed to the trans- 
formation block. This block consists of the quartz plate 
with fixed optical thickness S = 0.9046 and variable ori- 
entation a. So the state, which is to be measured, is de- 
termined by the parameter a. The measurement block is 
a Brown-Twiss scheme equipped with polarization filters 
placed in both arms (Fig. 3). Pulses coming from a cou- 
ple of single-photon modules (EG&G SPCM-AQR) were 
fed to the counter through a standard time-to-amplitude 
converter. 

In our experiments the exposure time for measuring 
each moment is 5 sec. This time is an important ex- 
perimental parameter. Each measurement consists of 
30 runs, after which the scheme is reset. Namely each 
measurement is performed by setting the angles of wave 
plates Xj 0j in both arms according to the tomo- 
graphic protocol (Table I). After 30 runs, a new set of 
angles is selected and the next moment is measured in 
the same way. The output data of the set-up are the 
mean coincidence counting rates. Examples of behavior 
for some moments {A, B, C, KcF, Im_F) versus the orien- 
tation of the plate QPl are plotted in Fig. 6. 



Measurement 



Protocol 1 




Transformation Preparation 

1 2,0) 



FIG. 5; Scheme of qutrits tomograph, consisting of three 
blocks. Preparation block includes pump laser(s) and non- 
linear crystal(s). Transformation block is the quartz plate 
(QPl) which orientation angle a determines the final state 
to be measured. The measurement block depends on the pro- 
tocol to be used (see Fig. 3 and Fig.4). 
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channel. Protocol 2 is easier to implement since only 
a single parameter n is changed whereas four Xi,2,^i,2 
parameters are varied in the Protocol 1. In perspective, 
this kind of protocol allows one to automate the quantum 
tomography procedure: the control plate can be rotated 
continuously and the reconstruction of the quantum state 
can be based on the analysis of coincidence rates corre- 
sponding to the respective values of = 1, .., 72). 



20 40 60 

Orientation of the quartz plate QPl, deg. 



D. Statistical reconstruction of biphoton- field 
qutrits from the outcomes of mutually 
complementary measurements 



Each of the 9 processes from Protocol 1 as well as of the 
72 processes from Protocol 2 is described by its ampli- 
tude M^. From the statistical point of view, the squared 
modulus of the process amplitude specifies the intensity 
of the event generation: 



FIG. 6: Some components of the matrix K4, versus the orien- 
tation of tlie quartz plate QPl. Different angles of the plate 
correspond to different states sent to the measurement block. 
The plot at the top corresponds to measured real moments 
A (squares), B (circles), C (triangles) and theoretical predic- 
tions A{—), B{ ), C{— o — ). The plot at the bottom shows 

measured complex moment KeF (squares), Im_F (triangles) 
and theoretical predictions ReF (-), Im_F (-). 



C. Experimental implementation. Protocol 2 

For the second method we used Ti:Sa laser with pulse 
duration about 250 fsec, operating at 800nm. After fre- 
quency doubling, the UV radiation with 400nm wave- 
length was sent into the same set-up as described above. 
For this protocol we used 2-mm BBO crystal cut for 
collinear degenerate type-I phase-matching. A quartz 
plate with the optical thickness S = 0.656 is placed after 
the BBO crystal to prepare the qutrit state to be mea- 
sured. The measurement part of the set-up was slightly 
changed (Fig. 4). Additional control quartz plate intro- 
duced in front of the beam splitter accomplishes the pro- 
tocol. Its orientation angle /i is a parameter defining 
the measurement process. The control plate is rotated 
with 5 degrees step from up to 360 degrees so that 
Protocol 2 consists of 72 measurements. Each arm of 
the Brown-Twiss scheme contains either quarter or half- 
wave plate with fixed orientation. The orientations are 
Xi = 18.8° for the quarter-wave plate in the first channel 
and 6*2 = —28.5° for the half-wave plate in the second 



= M*AU (22) 



The considered processes are examples of mutually com- 
plementary sets of measurements in the sense of Bohr's 
complementarity principle. The event-generation inten- 
sities Ri, for both protocols are the main quantities ac- 
cessible from the measurement. Making the bridge be- 
tween statistical and physical description of the process 
the quantities R^ coincide with the fourth moments in 
the field introduced above in Eq. (19). Their dimen- 
sion is frequency unit (Hz). The number of events oc- 
curring within any given time interval obeys the Poisson 
distribution. Therefore, the quantities i?^ specify the in- 
tensities of the corresponding mutually complementary 
Poisson processes and serve as estimates of the Poisson 
parameters (see below). 

Although the amplitudes of the processes cannot be 
measured directly, they are of the greatest interest as 
these quantities describe fundamental relationships in 
quantum physics. It follows from the superposition prin- 
ciple that the amplitudes are linearly related to the state- 
vector components. So the main purpose of quantum 
tomography is reproduction of the amplitudes and state 
vectors, which are hidden from direct observation. 

The linear transformation of the state vector c = 
{ci, 02,03} into the amplitude of the process M is de- 
scribed by a certain matrix X. For example, considering 
the first protocol this matrix can be easily obtained from 
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Table I (last column in Tabic I): 



X 



I l/^/2 





1/2 
1/2 
1/(2V2) 
V 1/(2V2) 




1/2 


l/(2v^) 
1/(2^2) 
-l/(2^/2) 
-*/(2V2) 





\ 


1/72 

-i/2 
-1/2 





«/(2%/2) 

-1/(2^/2) y 



(23) 



Using the matrix X the complete set of nine amplitudes 
of the processes can be expressed by a single equation, 



Xc^M 



(24) 



The matrix X is an instrumental matrix for a set of mutu- 
ally complementary measurements, by analogy with the 
conventional instrumental function. The implementation 
of the method to the first protocol has been considered 
in H. 

Consider an algorithm allowing one to calculate the in- 
strumental matrix X for Protocol 2. The matrix consists 
of 72 rows (the number of control plate orientations) and 
3 columns (the dimension of Hilbert space for qutrits). 
Each row is formed in the following way. Using coeffi- 
cients t^"'') and r("'') of the wave plates introduced to 
the signal and idler channels of the Brown-Twiss scheme 
(18a, b) the three-element row, which defines the process 
amplitude right after the control plate, can be written in 
the form 



A/4' A/2 



1 

71^ 



„(*) M 



fK^'''\'/A^\'/2 + '"a/2 ^ A/4 



''A/4''A/2 



(25) 

The unitary matrix G is defined by the control plate ac- 
cording to (7, 8), with a replacement a ^ /i, where /i is 
the control plate orientation (it takes 72 values from 0° to 
355° ). We chose the control plate to be a quarter- wave 
plate, so S = 7r/4. Finally, 



G = G(mO,*-1,2, 



,72. 



(26) 



Each row of the instrumental matrix X (72 rows, 3 
columns) is defined by the product of the row I (which 
is the same for any process) and the matrix G (which is 
defined by the control plate orientation angle) 



X,=lG{fi,),i = 1,2,. ..,72, 
where Xi is the ith row of the matrix X . 



(27) 



of experimental data is limited in this case, the recon- 
structed density matrix may have non-physical proper- 
ties like negative eigenvalues. But in the general case 
of s— dimensional systems the problem of density matrix 
reconstruction using direct results of measurements can 
not be solved since the corresponding inverse problem is 
ill-posed. 

When analyzing the experimental data, we use the so- 
called root estimator of quantum states [2^. This ap- 
proach is designed specially for the analysis of mutu- 
ally complementary measurements (in the sense of Bohr's 
complementarity principle). The advantage of this ap- 
proach consists of the possibility of reconstructing states 
in a high- dimensional Hilbert space and reaching the ac- 
curacy of reconstruction of an unknown quantum state 
close to its fundamental limit. Below we consider two 
methods of quantum state root estimation that give sim- 
ilar results. They are the least-squares method (LSM) 
and maximum- likelihood method (MLM). 



A. Least-squares method 

In statistical terms, Eq. (24) is a linear regression 
equation. A distinctive feature of the problem is that 
only the absolute value of the process amplitude M is 
measured in the experiment. The estimate of the abso- 
lute value of the amplitude is given by the square root of 
the corresponding experimentally measured coincidence 
rate: 



(28) 



where k^, is the number of events (coincidences) detected 
in the I'-th process during the measurement time t. 

It is worth noting that, by the action of the root- 
square procedure on a Poissonian random value, one gets 
at the random variable with a uniform variance, i.e., at 
the variance stabilization [23|. Note also, since we deal 
not with event probabilities but with their rates or in- 
tensities, it is convenient to use un-normalized state vec- 
tors. These vectors allow the coincidence counting rate 
(event-generation intensities) to be derived directly from 
Eqs. (22), (24) without introducing the coefficients re- 
lated to the biphoton generation rate, detector efficien- 
cies, etc. The dimensionality of the vector state obtained 
in such a way is 1/ vtimc. The final state vector obtained 
by the reconstruction procedure, nevertheless, should be 
normalized to unity. 

Assuming that the variances of different |Afi,|™P are 
independent and identical, one can apply the standard 
least- squares estimate to Eq. (24) 



I app 

m 



IV. METHODS OF QUANTUM STATE 
RECONSTRUCTION 

In the simplest case the density matrix can be esti- 
mated directly from the measurements. Since the set 



c= {X^'Xy^X^'M (29) 

Unlike the traditional least-squares method, Eq. (29) 
cannot be used for the explicit estimation of the state vec- 
tor c, because it is to be solved by the iteration method. 
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The absolute value of M is known from the experiment 
(I My I = |Mi,|™P). We assume that the phase of vector 
Xc at the i-th iteration step determines the phase of the 
vector M at the i + 1-th step. In other words the phase 
is determined by the iteration procedure. 

It turns out that, for the Gaussian approximation of 
Poisson's quantities, this least-squares estimate coincides 
with a more exact and rigorous maximum-likelihood es- 
timate considered below. 



B. Maximum-likelihood method 

The likelihood function is defined by the product of 
Poissonian probabilities: 

i 

where ki is the number of coincidences observed in the 
i-th process during the measurement time t^, and are 
the unknown theoretical event-generation intensities (ex- 
pected number of coincidences proportional to the mo- 
ments in the field), whose estimation is the subject of 
this section. 

The logarithm of the likelihood function is, if we omit 
an insignificant constant, 

InL = ^ [h HKU) - KU). (31) 

i 

Let us introduce the matrices with the elements defined 
by the following formulas: 

Ija ^'^^tiX*jXis, (32) 



•^j«=Ev^u-^-' .7-s = l>2,3. (33) 

The matrix / is determined from the experimental proto- 
col and, thus, is known a 'priori (before the experiment). 
This is the Hermitian matrix of Fisher's information. The 
matrix J is determined by the experimental values of ki 
and by the unknown event-generation intensities . This 
is the empirical matrix of Fisher's information (see also 
Appendix) . 

In terms of these matrices, the condition for the ex- 
trcmum of the function (31) can be written as 

Ic ^ Jc. (34) 

Hence, it follows that 

/-I Jc = c. (35) 



The latest relationship is known as the likelihood equa- 
tion. This is a nonlinear equation, because Xi depends 
on the unknown state vector c. Because of the simple 
quasi-linear structure, this equation can easily be , solved 
by the iteration method j26| . The quasi- identity opera- 
tor /"^ J acts as the identical operator upon only a single 
vector in the Hilbert space, namely, on the vector corre- 
sponding to the solution of Eq. (35) and representing the 
maximum possible likelihood estimate for the state vec- 
tor. The condition for the existence of the matrix is 
a condition imposed on the initial experimental protocol. 
The resulting set of equations automatically includes the 
normalization condition, which is written as 

i i 

This condition implies that, for all processes, the total 
number of detected events is equal to the sum of the 
products of event detection frequencies during the mea- 
surement time. 



C. Analysis of the experimental data 

1. Pure state reconstruction 

The examples of qutrit state reconstruction using both 
the least-squares and maximum-likelihood methods are 
given in Table II. 

The value of the fidelity parameter F is defined as 

F — \ {^theory \^exp)\ ■ {^^) 

It gives the conventional measure of correspondence be- 
tween the theoretical and experimental state vectors. 

The dependence of fidelity on the amount of exper- 
imental data obtained is shown in Fig. 7. This figure 
shows the fidelity achieved in the experiment in compar- 
ison with the theoretical range (see Appendix for more 
details). The lower boundary corresponds to 5% quan- 
tile of statistical distribution, while the upper to 95%- 
quantile. It is clearly seen that the fidelity value achieved 
experimentally for a small volume of experimental data 
is completely within the limits of the theoretical range, 
while it goes out for a higher volume. Such behavior of 
fidelity is due to the existence of two different error types 
arising under the reconstruction of quantum states. Let 
us call them statistical and instrumental errors, respec- 
tively. The statistical errors are caused by a finite number 
of quantum system representatives to be measured. As 
the measurement time increases, the information about 
the quantum state of interest progressively increases (see 
Appendix). Accordingly, the statistical error becomes 
smaller. The instrumental errors are caused by the re- 
searcher's incomplete knowledge of the system; i.e., more 
exact information exists, in principle, but it is inacces- 
sible to the experimenter. Thus, a comparison between 
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TABLE II: Results of the state reconstruction. The left column indicates the orientation of the quartz plate QPl, determining 
the state to be measured. Values of the optical thickness of QPl are 5 = 0.656 for the pulsed regime (Protocols 1, 2), and 
S = 0.9046 for the cw regime (Protocol 1). Theoretical state vectors are placed in the right column. The table contains the 
amplitudes of the reconstructed states (ci,C2,C3) as well as their fidelities, calculated by least-squared (LSM) and maximum- 
likehhood (MLM) methods. 



Pulsed regime, 5 — 0.656, Protocol 1 


a 


Fids 
LSM 


5lity 
MLM 


State vector 
LSM 


experiment 
MLM 


State vector: theory 

{ci,C2,C3)theory 


0° 


0.99981 


0.99979 


-0.0046+0.00401 
-0.0050-0.01151 
0.9999 


-0.0065+0.00571 
-0.0053-0.01021 
0.9999 





1 


40° 


0.9989 


0.9989 


-0.3669-0.06911 
-0.0657-f0.6814i 
0.6261 


-0.3669-0.06871 
-0.0653+0.68151 
0.6261 


-0.3482-0.09481 
-0.0900+0.67321 
0.6392 


80° 


0.9993 


0.9993 


-0.0088-f0.0439i 
0.16914-0.25871 
0.9500 


-0.0091+0.04391 
0.1697+0.2591 
0.9498 


-0.0136+0.04131 
0.1691+0.23381 
0.9565 


Pulsed regime, 5 — 0.656, Protocol 2 


0° 


0.99846 


0.99847 


-0.0071-0.01351 
0.03594-0.00461 
0.9992 


-0.0072-0.01351 
0.0357+0.00461 
0.9992 




1 


40° 


0.9991 


0.9991 


-0.3442-0.11391 
-0.0987-f0.6546i 
0.6560 


-0.3444-0.11421 
-0.0990+0.65451 
0.6559 


-0.3482-0.09481 
-0.0900+0.67321 
0.6392 


80° 


0.9981 


0.9981 


-0.0093+0.04301 
0.2122+0.24081 
0.9461 


-0.0094+0.04301 
0.2121+0.24081 
0.9461 


-0.0136+0.04131 
0.1691+0.23381 
0.9565 


CW-regime, S = 0.9046, Protocol 1 


0° 


0.99325 


0.99313 


-0.0030-0.05121 

0.9966 
-0.0015-0.06421 


-0.0028-0.05141 

0.9966 
-0.0013-0.06491 



1 



60° 


0.9886 


0.9799 


0.7236 
0.1165-0.12311 
0.2792+0.60801 


0.7244 
0.1245-0.12101 
0.1694+0.64531 


0.7052 
0.0392-0.06161 
0.2990+0.63871 



Thus, for a small volume of experimental data, statis- 
tical errors prevail, whereas for large sample sizes, the 
setting errors and the instability of protocol parameters 
dominate. The number of events at which the statistical 
error becomes smaller than the instrumental error can 
be called the coherence volume. Numerically coherence 
volume can be estimated as the intersection point be- 
tween the experimental fidelity and the lower theoretical 
fidelity curve. In our case this value is about 25,000- 
30,000 events. Starting approximately from this value, 
fidelity is reaching saturation and further growth of ex- 
perimental data volume does not lead to an increase in 
the precision of quantum system estimation. 

Fig. 7 relates to the state defined by orientation angle 
of quartz plate QPl a ^ 50° (for Protocol 2). To plot 
Fig. 7 we used the following technique for passing from 
full- volume experiment to a partial- volume experiment. 
Let us consider the parameter < / < 1 that charac- 
terizes the volume of experimental data. Suppose that 
/ = 1 for a full- volume experiment. A partial volume 
experiment may be introduced considering the observa- 
tion time t'j^ — ftu instead of Hence, performing 
a single full-volume experiment means providing with a 
large (practically infinite) number of partial-volume ex- 




1000 10000 100000 1000000 

Number of events 

FIG. 7: Fidelity dependence on the sample size. Mean values 
and standard deviations corresponding to the sample volumes 
/ = 0.01; 0.04; 0.1; 0.25; 0.5; 1.0. 



the state reconstruction result and the fundamental sta- 
tistical level of accuracy can serve as a guide for the pa- 
rameter adjustment of the set-up. 
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FIG. 8: Informational x'^ criterion for small sample 
sizes:sample size=4400. 



periments. 

For a given volume of experimental data / each event 
from the full-volume experiment is picked up with the 
probability / and rejected with the probability 1 — f. 
Due to the presence of statistical fluctuations the equa- 
tion for the number of observations. k^{t'j^) = fki,(tu), is 
violated. Therefore a unique estimate of the state vector 
corresponds to every partial-volume experiment. Fig. 7 
shows mean values and standard deviations correspond- 
ing to volumes / = 0.01; 0.04; 0.1; 0.25; 0.5; 1.0. For 
each / < 1 ten experiments were simulated. 

The results of informational fidelity research, intro- 
duced in Appendix, are shown in Figs. 8, 9. These figures 
correspond to the same data set as shown in Fig. 7. Distri- 
bution density of informational fidelity for a small (com- 
pared to the coherence volume) sample size closely agrees 
with the theoretical result given by (A14) (see Fig. 8). In 
this case the instrumental error is negligibly small com- 
pared to the statistical one. When the sample volume is 
close to the coherence volume (Fig. 9) the influence of in- 
strumental and statistical errors is about equal. In other 
words, the informational losses caused by averaging over 
instrumental errors are approximately equal to the losses 
caused by statistical ones. Finally, if the sample size is 
greater than the coherence volume, instrumental errors 
predominate. It means that the statistical informational 
errors arc negligibly small compared to the instrumental 
ones. 



2. Mixture separation algorithm 



FIG. 9: Informational criterion for large sample sizes: sam- 
ple size=27750. The disagreement between observations and 
theoretical curve for large sample sizes is due to the instru- 
mental error. 



intensity. 



(2) 



(38) 



where u — 1,2,..., Vmax and I'max is the total number of 
processes, A^^'' and a!?'' are the estimates of intensities of 
processes for a given step of the iteration procedure. 

At a certain iteration step, let us represent k,^ as a sum 
of two components. 



(39) 



For each component, we can obtain the estimates for the 
state vector, amplitudes, and intensity of the processes 
according to the method of pure state analysis described 
in the previous section. Since we get new intensity es- 
timates, let us again split the total number of events in 
every process proportionally to the intensities of the com- 
ponents. In such a way, a new iteration is formed and 
the whole procedure is repeated. The described process 
is called quasi-Bayesian algorithm p^ . 

As a result, the iteration process converges to some 
(non-normalized) components c*-^' and c^^\ Thus, the 
mixture separation algorithm reduces to numerous es- 
timations of pure components according to the simple 
algorithm described above in section IV. B. As a result 
of the whole algorithm execution, the estimate for the 
density matrix of the mixture appears: 



Let us describe the algorithm for reconstructing a two- 
component mixed state. This algorithm can be easily 
generalized to an arbitrary number of components. 

The total number of events observed in every process 
is divided between the components proportional to the 



p = cWc(i)t+c(2)c(2)t, 



Tr(p)- 



(40) 



(41) 
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TABLE III: Example of the mixture separation using quasi- 
Bayesian algorithm for the given state. CW-regime, protocol 
1. 



State vector:theory 


Density matrix: experiment 


a = 30° 
5 = 0.9046 

(Cl , C2 5 C3 ^theory 


First principal 

component 
weight=0.9238 

(Ci , C2 , C3 ) exp 


Second principal 

component 
weight=0.0762 


0.7052 
-0.0392-0.06161 
0.2990-0.63871 


0.7019 
-0.0466-0.13251 
0.2245-0.66121 


-0.3027-0.28581 
-0.6529-f0.3291i 
0.5140-0.16681 


Fidelity=0.9916 



The last procedure is normalization of the density matrix. 

A remarkable feature of the algorithm is that ac- 
cording to numerical calculations, independent on zero- 
approximation selection of the mixture components, the 
resulting density matrix p is always the same. Of course, 
the components c^^-* and c^^-* are different for the random 
selection of the zero-approximation. 

The mixed state reconstruction accuracy is described 
by the following fidelity: 



F 



(42) 



where ' and p are the exact and reconstructed 
density matrices respectively. For a pure state 
(p2 ^ (p(0))2 ^ ^(0)) fidelity (42) converts to (37). 

Actually in the present work we did not intend to 
generate a given mixed state of qutrit in experiment. 



it will be done later Nevertheless, applying the 

described algorithm to the data we can check whether 
the state produced in our system is pure. For example, 
consider the case when the state ^2 = is fed to 

the quartz plate QPl (see Table II). This state is the 
most interesting to be tested, since \H,V) and \V,H) 
are distinguishable due to the polarization dispersion in 
BBO crystal. Namely, extraordinarily polarized photons 
{H) propagate faster than ordinary {V) ones in the crys- 
tal. Therefore a group velocity compensator has to be 
used for making them indistinguishable |23 . Non-perfect 
compensation (we reached 95% visibility for polarization 
interference) is the main reason why the fidelity recon- 
struction for these states is not so high. The results of 
applying quasi-Bayesian algorithm to the reconstructed 
state arc in Table III. Wc chose the state correspond- 
ing to the angle a — 30°. It is clearly seen that the 
weight of the first principal component is much greater 
that of the second one. Doing the same procedure with 
the ^I*! = |2, 0) initial state, we have checked that the es- 
timator for a pure state vector is extremely close to the 
estimator of the major density matrix component. 

To illustrate the quasi-Bayesian approach, let us con- 
sider a result of reconstruction for a two-component mix- 
ture using protocol 1. Suppose one has a mixture of two 
pure states prepared from |2,0) by quartz plates QPl' 
and QPl" oriented at angles a = —30° and a = 50°, 
respectively. Let the optical thickness of both plates be 
5 = 0.656. Ten thousands events were generated (on the 
average) for every component. The theoretical density 
matrix for the mixed state under consideration is 



/ 0.1134 0.0263 -f0.0808i -0.1987 - 0.0558i 

= 0.0263 - 0.0808i 0.4404 0.0679 + 0.0752i | (43) 

\ -0.1987 + 0.0558i 0.0679 - 0.0752i 0.4462 

A typical example of a reconstructed density matrix is the following: 



0.1162 
0.0204- 0.0808i 
-0.1965 -H0.0691i 



0.0294 + 0.0808i 

0.4298 
0.0697- 0.0796i 



-0.1965 - 0.069H' 
0.0697 + 0.0796i 
0.4540 



(44) 



The reconstructed matrix fidelity is = 0.999431. Anal- 
ysis of the principal components of density matrix is given 
in Table IV. 

This example shows a reasonably high accuracy of 
mixed state reconstruction. The statistical properties of 
the proposed algorithm were studied by means of the 
Monte-Carlo method. One hundred numerical experi- 
ments were conducted similar to the one described above. 
To verify the reliability, the solution was found twice for 
each experiment (with random zero approximation selec- 



tion). The solutions appeared to be equal for all cases 
(within negligible small computational error). The ob- 
tained statistical fidelity distribution is shown in Fig. 10. 
Numerical research shows that the fidelity distribution 
density is well described by the beta-distribution. 
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TABLE IV: Analysis of the principal components of the den- 
sity matrix for the state (43, 44): numerical simulation. 



State vector 

(Cl, C2, C3) 


Fidelity 


First principal component 


Experiment 
weight=0.6188 


Theory 
weight=0.6143 




-0.3658-0.04481 
0.2085-^0.47431 
0.7718 


-0.3668-0.02111 
0.2294+0.49341 
0.7543 


0.9985 


Second principal component 


Experiment 
weight=0.3812 


Theory 
weight=0.3857 




-0.1208-0.26431 
-0.1659-0.81501 
0.4731 


-0.1490-0.23821 
-0.1986-0.79421 
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FIG. 10: Simulation of the fidelity between theoretical and 
reconstructed density matrices in a mixture separation prob- 
lem. 100 numerical experiments of 20,000 events per each (on 
the average) were made. 

V. CONCLUSION 

The procedure of quantum state measurement for a 
three-state optical system formed by a frequency- and 
spatially degenerate two-photon field has been consid- 
ered in this work. The method of the statistical estima- 
tion of the quantum state through solving the Ukelihood 
equation and examining the statistical properties of the 
resulting estimates has been developed. Based on the ex- 
perimental data (fourth-order moments in the field) and 
the root method of estimating quantum states, the initial 
wave function of qutrits has been reconstructed. 

Experimental data analysis is based on representing 
the event generation intensity for each one of mutually- 
complementary quantum processes as a squared module 
of some amplitude. A complete set of measured pro- 
cesses amplitudes can be compactly described using the 
instrumental matrix. In the framework of the formal- 
ism of a process amplitude one can apply effective tools 
for the quantum state reconstruction: least-squares and 



maximum-likelihood methods. 

The developed analysis tools provide means of quan- 
tum state reconstruction from the experimental data 
with high accuracy and reliability. The estimate accu- 
racy is determined by concurrence of two types of errors: 
statistical ones and instrumental ones. For smaller sam- 
ple sizes statistical errors are dominant, while for greater 
ones instrumental errors dominate. 

Instrumental errors lead to fidelity saturation at less 
than unity level. In the present work, fidelity for most of 
performed experiments (more than 20) exceeded level of 
0.995. For many cases the level of 0.9998 was achieved. 
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APPENDIX A: STATISTICAL FLUCTUATIONS 
OF THE STATE VECTOR 

As it was already mentioned above, an un-normalizcd 
state vector provides the most complete information 
about a quantum system. The usage of un-normalizcd 
vector allows us to remove an interaction constant in Eq. 
(24). The norm of the vector c, obtained as a result of 
quantum system reconstruction, provides one with the 
information about the total intensity of all the processes 
considered in the experiment. However, the fluctuations 
of the quantum state (and norm fluctuations, in particu- 
lar) in a normally functioning quantum information sys- 
tem should be within certain range defined by the sta- 
tistical theory. The present section is devoted to this 
problem. 

Practical significance of accounting for statistical fiuc- 
tuations in a quantum system relates to developing meth- 
ods of estimation and control of precision and stability 
of a quantum information system evolution, as well as 
methods of detecting external interception (Eve's attack 
on the quantum channel between Alice and Bob). 

The estimate of the un-normalized state vector c, ob- 
tained by the maximum-likelihood principle, differs from 
the exact state vector c^"' by random values 5c = c'-''^ — c. 
Let us consider the statistical properties of the fluctua- 
tion vector 5c by expansion of the log likelihood function 
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near the stationary point, 



SlnL 



(Al) 

Together with the Hermitian matrix of the Fisher infor- 
mation / (32), we define the symmetric Fisher informa- 
tion matrix K, whose elements are defined by the foUow- 
ing equation: 



(A2) 



where M^, is the amphtude of the v- th process. In the 
general case, K is a complex symmetric non-Hermitian 
matrix. From all possible types of fluctuations, let us pick 
out the so-called gauge fluctuations. Infinitesimal global 
gauge transformations of a state vector are as follows: 



iecjj = 1,2, ...,s 



(A3) 



where e is an arbitrary small real number, s is the Hilbert 
space dimension. 

Evidently, for gauge transformations SlnL = 0. It 
means that two state vectors that differ by a gauge trans- 
formation, are statistically equivalent, i.e. they have 
the same likelihood. Such vectors are physically equiv- 
alent since the global phase of the state vector is non- 
observable. From statistical point of view, the set of mu- 
tually complementing measurements should be chosen in 
such a way that for all other fluctuations (except gauge 
fluctuations) dliiL < 0. This inequality serves as the 
statistical completeness condition for the set of mutually 
complementing measurements. 

Let us derive some constructive criteria of the statisti- 
cal completeness of measurements. The complex fluctua- 
tion vector Sc is conveniently represented by a real vector 
of double length. After extracting the real and the imag- 
inary parts of the fluctuation vector Scj ~ Jc^^^ + iScf^ 
we transfer from the complex vector Sc to the real one 
S^: 



fSc['^\ 



Sc 



Sci 



SCs 



S^ = 



Sc 
Sc 



(1) 

s 

(2) 



(A4) 



In the particular case of qutrits (s = 3) this transition 
provides us with a 6-component real vector instead of a 
3-component complex vector. 

In the new representation Eq. (Al) becomes: 



SlnL ^ -Hs.St^S^j ^ -{S(\H\SO, 



(A5) 



where matrix H is the "complete information matrix" 
possessing the following block form: 



f Re(/ + K) -Im(/ + K) 
\ Im(/ - K) Re(/ - K) 



(A6) 



The matrix H is real and symmetric. It is of double 
dimension respectively to the matrices / and K. For 
qutrits, / and K are 3x3 matrices, while iJ is 6 x 6. 

Using matrix H it is easy to formulate the desired char- 
acteristic completeness condition for a mutually comple- 
menting set of measurements. For a set of measurements 
to be statistically complete, it is necessary and sufficient 
that one and only one eigenvalue of the complete infor- 
mation matrix H is equal to zero, while the other ones 
are strictly positive. 

We would like to stress that checking the condition one 
not only verifies the statistical completeness of a mea- 
surement protocol but also, insures that the obtained ex- 
tremum is of maximum likelihood. 

An eigenvector that has eigenvalue equal to zero corre- 
sponds to gauge fiuctuation direction. Such fluctuations 
do not have physical meaning as stated above. Eigenvec- 
tors corresponding to the other eigenvalues specify the 
direction of fluctuations in the Hilbert space. 

The principal fluctuation variance is 



1 

2^' 



l,...2s-l, 



(A7) 



where hj is the eigenvalue of the information matrix H, 
corresponding to the j-th principal direction. 

The most critical direction in the Hilbert space is 
the one with the maximum variance cr|, while the cor- 
responding eigenvalue hj is accordingly minimal. The 
knowledge of the numerical dependence of statistical fluc- 
tuations allows one to estimate distributions of various 
statistical characteristics. 

The important information criterion that specifies the 
general possible level of statistical fiuctuations in quan- 
tum information system is the chi-square criterion. It 
can be expressed as 



2(S^\H\SO'^X^{2s-l), 



(A8) 



where s is the Hilbert space dimension 

The left-hand side of Eq. (A8), which describes the 
level of state vector information fluctuations, is a chi- 
square distribution with 2s — 1 degrees of freedom. 

Validity of the analytical expression (A8) is justified 
by the results of numerical modeling and observed data. 
Similarly to Eq. (A4), let us introduce the transforma- 
tion of a complex state vector to a real vector of double 
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length: 



The mean value of relative information fluctuations is 



Cl 



J2) 



(A9) 



It can be shown that the information carried by a state 
vector is equal to the doubled total number of observa- 
tions in all processes: 



m\0 = 2r 



(AlO) 



where n = Y^k^. 



Then, the chi-square criterion can be expressed in the 
form invariant to the state vector scale (let us remind 
that we consider a non- normalized state vector). 



m\i) 



An 



(All) 



Relation (All) describes the distribution of relative in- 
formation fluctuations. It shows that relative informa- 
tion uncertainty of a quantum state decreases with the 
number of observations as 1/n. 



2s -1 
An 



(A12) 



The information fidelity may be introduced as a mea- 
sure of correspondence between the theoretical state vec- 
tor and its estimate: 



H 



{mm 



(A13) 



Correspondingly, the value I^Fh is the information loss. 

The convenience of Fh relies on its simpler statistical 
properties compared to the conventional fidelity F. For 
a system where statistical fiuctuations dominate, fidelity 
is a random value, based on the chi-square distribution. 



H 



1 



An 



(AM) 



where x^(2s — 1) is a random value of chi-square type 
with 2s — 1 degrees of freedom. 

Information fidelity asymptotically tends to unity 
when the sample size is growing up. Complementary 
to statistical fluctuations noise leads to a decrease in the 
informational fidelity level compared to the theoretical 
level (A14). 
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